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We construct, for the first time, a model of graceful exit transition from a dilaton-driven infla- 
tionary phase to a decelerated Friedman— Robertson— Walker era. Exploiting a demonstration that 
classical corrections can stabilize a high curvature string phase while the evolution is still in the 
weakly coupled regime, we show that if additional terms of the type that may result from quantum 
corrections to the string effective action exist, and induce violation of the null energy condition, then 
evolution towards a decelerated Friedman— Robertson— Walker phase is possible. We also observe 
that stabilizing the dilaton at a fixed value, either by capture in a potential minimum or by radiation 
production, may require that these quantum corrections are turned off, perhaps by non-perturbative 
, effects or higher order contributions which overturn the null energy condition violation. 
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I. INTRODUCTION 



< 
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An inflationary scenario , inspired by duality symmetries of string cosmology equations , is based on the 

^ fact that cosmological solutions to string dilaton-gravity come in duality-related pairs, the plus branch and the 
\^ minus branch (— ) The (-I-) branch has kinetic inflationary solutions in which the Hubble parameter increases 
I with time. The minus branch (— ) can be connected smoothly to a standard Friedman— Robertson— Walker (FRW) 
■ decelerated expansion of the Universe with constant dilaton. The scenario (the so called "pre-big-bang" scenario) is 
' that evolution of the Universe starts from a state of very small curvature and coupling and then undergoes a long 
phase of dilaton-driven kinetic inflation described by the {+) branch and at some later time joins smoothly standard 
radiation dominated cosmological evolution, thus giving rise to a singularity free inflationary cosmology. Recently, 
the required initial conditions of the Universe in this scenario were discussed . Our focus in this paper is mainly 
i^H ■ on the flnal fate of the Universe at late times and we simply assume that the correct initial conditions were chosen 
I ' such that a long dilaton-driven inflationary phase was indeed part of the evolution. 
Oh! 

The graceful exit transition from the initial long phase of dilaton-driven kinetic inflation to the subsequent standard 
' radiation dominated evolution has been a subject of many investigations [[7|-p^. In |^ it was argued, and later proved 
• • , PI , that such a transition cannot occur while curvature was below the string scale and the string coupling was still 
. ^ weak, leading to the conclusion that an intermediate "string phase" of high curvature (previously suggested as a 
possibility ||l|Q) or strong coupling is actually required fl^ . 
^ ' In we proposed to use an effective description in terms of sources that represent arbitrary corrections to the 
[ lowest order equations. We were able to relate necessary conditions for graceful exit to energy conditions appearing 
in singularity theorems of Einstein's general relativity | p^ . In particular, we showed that a successful exit requires 
violations of the null energy condition (NEC). Since most classical sources obey NEC this conclusion hints that 
quantum effects, known to violate NEC in some cases, may be the correct sources to look at. To briefly recap the 
relevant results of we need to recall that there are different conformal frames in which to describe the equations. 
These are related by local field redefinitions which, supposedly, do not affect physical observables The two frames 
that we use are the string frame (5J_ ) and the "lowest order Einstein frame" (_E[ ) in which the lowest order classical 
action kinetic terms of gravity and the dilaton are diagonal. For brevity we call _E[ the Einstein frame. The general 
analysis of |^] resulted in a set of necessary conditions on the evolution in terms of the Hubble parameters Hs in the 
string frame and He in the Einstein frame and the dilaton 0. We include them for completeness 

• Initial conditions of a (+) branch and Hs, (j) > require He < 0. 

• A branch change from (+) to (— ) has to occur while He < 0. 

• A successful escape and exit completion requires NEC violation accompanied by a bounce in El after the branch 
change has occurred, ending up with He > 0. 

• Further evolution is required to bring about a radiation dominated era in which the dilaton effectively decouples 
from the "matter" sources. 
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The question we set to answer in this paper is the following. Suppose that effective sources of the type that are 
expected to appear as corrections to the lowest order effective action of strings do provide NEC violation. Then would 
a complete exit transition actually occur? We answer this question in the affirmative, but some surprising obstacles 
are found on the way and we also suggest ways to overcome these obstacles. 

Because in our scenario the Universe evolves towards higher curvatures and stronger coupling, there will be some 
time when the lowest order effective action can no longer reliably describe the dynamics and it must be corrected. 
Corrections to the lowest order effective action come from two sources. The first are classical corrections, due to 
the finite size of strings, arising when the fields are varying over the string length scale As = \fo! . These terms 
are important in the regime of large curvature. The second are quantum loop corrections. The loop expansion is 
parameterized by powers of the string coupling parameter e'^ = g^tringi which in the models that we consider is, of 
course, time dependent. So quantum corrections will become important when the dilaton becomes large, the regime 
we refer to as strong coupling. 

In this paper we make use of both types of corrections. The role of a' corrections [^| is to create an attractive fixed 
point which stabilizes the evolution in a high curvature regime with linearly growing dilaton. The basin of attraction 
of this fixed point is large enough to allow it to be reached from generic initial conditions. Further, the location of 
this fixed point is such it forces the evolution to undergo a branch change, and all of this may occur for small values 
of the dilaton (weak coupling) so the quantum corrections can be ignored. With the linearly growing dilaton, the 
quantum corrections will eventually become important. It is these we will attempt to use to let the Universe escape 
the fixed point and complete the transition to a decelerated FRW evolution. We allow ourselves the freedom to choose 
the coefficients of the quantum correction terms arbitrarily, in particular, their sign is chosen so as to induce NEC 
violation. Our reasoning for allowing this freedom stems in part from a lack of any real string calculations and in part 
by our desire to verify by constructing explicit examples the general arguments of 0. 

We use the general framework set up in, and rely for analytical calculations mostly on We perform numerical 
integration of specific equations with limited new analytical considerations. We show, by constructing explicit exam- 
ples, that a complete exit transition is not forbidden, which, of course, does not yet prove that in string theory it 
actually occurs. For that, an explicit string computation of quantum correction terms is necessary. This may also 
require taking into account the backreaction of the particles produced during the dilaton-driven phase. 

A by-product of our analysis is a detailed description of the high-curvature ( "string" ) phase in between the dilaton- 
driven inflationary phase and the decelerated FRW phase. We obtain a more complicated phase than the one postu- 
lated in m with non-constant </> and Us- We do not take the details of the evolution in the string phase too seriously, 
because the terms that determine those details have arbitrary coefficients. However, our examples could be taken as 
an indication of what the real string phase may eventually look like. 

The paper is organized as follows. In section II we present and discuss solutions which exhibit branch change and 
exit completion with a' and quantum corrections. In section III we explain the remaining problem of correction 
dominated evolution, show how to model a shut-off of the corrections and verify that once corrections are shut off, 
the usual decelerated evolution such as radiation domination or capturing of the dilaton in a minimum of a potential 
may follow. We summarize our results and explain which further calculations are necessary to reach a conclusive 
statement about exit transition in string models. In a technical appendix we outline details of derivation of equations 
of motion, numerical integration etc. 



II. A MODEL FOR EXIT COMPLETION 



String theory effective action in four dimensions takes the following general form. 
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(2.1) 



where g^jy is the 4-d metric and is the dilaton, the effective action is written here in the string frame (5J. ). The 
Lagrangian Lc may contain arbitrary corrections to the lowest order 4-d action coming from a variety of sources as 
specified below. 



We are interested in solutions to the equations of motion derived from the action (2.1) of the FRW type with 
vanishing spatial curvature (non-vanishing spatial curvature may be included as an effective source in the equations) 
ds^ = —dt\ -I- Ogii^dxidx^ and — <j>{t). We will assemble the equations of motion by deriving the correction energy- 
momentum tensor T^^ = ^^^^fa, which will have the form T'^ = diag{p, —p, —p, —p). In addition we have 

1 Sy/^Cc 



another form of source term arising from the variation by 
In terms of these sources the equations of motion are 



equation, A^Cc 
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ml + - iHs^ = ie*p (2.2) 



-2Hs - iHi, + 2i7s<?!' - \r + <j) = ie> (2.3) 
3ijs + 6i/|-3i/s0+i(/)2_0= le^A^^e (2.4) 
p + iHs{p + p)^-/^^Cc^, (2.5) 



= o-s/oLSi E^nd we have fixed our units such that IGna' = 1. 
We will make the split 

Cc ^ Ca' + Cg + (2.6) 

and we will divide up the source terms analogously. For example, p = p^' + Pq + Pm- We take, for concreteness, the 
Lagrangian to be of the specific form of the a' classical corrections proposed in [|l^ , 

K«-e-n%-(^). (2.7) 

The contribution of this term to the effective sources is 

p^,=e-^{6Hl^~^) (2.8) 

Pa' = e-'^(-4i?|0 - iHsHscj) + 2i7|02 2F|,/,) (2.9) 

3(i** ■ •• 

A^/:„, = e-'^i-GHl - mlHs + 3Hs(f>^ " ^ + (2-10) 

Although we have fixed our units differently from [l0|] (where ka' = 1, k depending on the species of string theory), 
we have chosen the coefficient of this term so as to lead to solutions which arc numerically identical to theirs for 
purposes of comparison. 

The Lagrangian Cq will designate quantum loop corrections parameterized by powers of e'^. We will examine a 
variety of candidates for Cq. The Lagrangian £,„ will represent such things as radiation or a dilaton potential when 
we discuss the final transition to radiation domination and a fixed dilaton. 

Setting Cq ~ £,„ = 0, for the moment, we integrate the equations of motion starting from initial condit ions near 



the (+) branch vacuum. To do this we solve (2_^) and {2A) for t he h ighest derivatives of Hs and (f> and use (2.2) as a 



constraint on the initial conditions. The conservation equation (2.5) is an identity owing to the usual redundancy in 
the equations, only when we add radiation to the system will we use this equation in the evolution. We recover the 
results of in four dimensions where they find attraction into a fixed point at Hs = 0.616..., </> = 1.404.... 

As shown in Fig. |l|a the solution begins near the (+) branch vacuum and makes a branch change but does not 
"complete" the exit by proceeding to the p > region. The quantity characterizing branch sign is plotted in Fig. |l|e. 
We examined more general forms of this Lagrangian involving different coefficients. Varying the coefficients moves 
the fixed point. For some ratios of coefficients the solutions encounter singularities before reaching the fixed point 
(we discuss the form of these singularities in section 3) but other ratios yield solutions of the same generic type as we 
have presented here. 

In spite of considerable progress we still need to look more closely at the sources and Einstein frame evolution to 
answer the question of how close we have come to the true goal of decelerated FRW evolution. 

The evolution of p and p is easily found by substituting the solution back into the equations of motion and solving 
for the sources. We see that this phase does violate NEC in the string frame. However, recalling our stronger claim 
that there should also be NEC violation in the "lowest order Einstein frame" (denoted with a subscript E and as 
explained in the introduction, many times for brevity called just the Einstein frame), we can compute the value of 
the Einstein frame sources by finding the Hubble parameter He and its derivative He, 

He = e^/^Hs ~ ^<P) (2.11) 

He = e^{Hs - + i0i?s - ^) (2.12) 

All quantities on the right hand side of these are given strictly in terms of string frame quantities and time. The 
Einstein equations are, 
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Notice that we have absorbed the dilaton kinetic energy into the definition of pE and pE- 



(2.13) 
(2.14) 
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FIG. 1. A solution to the classical equations with only a' corrections, Cq 



(f) 

0. (a) Evolution in the (0/3, i?s) plane. 



The four lines plotted are, in order of increasing slope, the (+) branch vacuum (p = 0, <?!> = (3 + 3^^^)Hs), branch change 
((/) = 3Hs), Einstein frame "bounce" {(p = 2Hs) and ( — ) branch vacuum (p = 0,</>= (3-3'''2)ii's). The cross is at the location 
of the fixed point. The remaining figures are various quantities plotted as a function of string time, (b) The Einstein frame 
Hubble parameter He- (c) The string frame source terms e'^{ps +Ps) and e'^ps- (d) The Einstein frame source term pE +pe- 
(e) (f> — SHs, a quantity indicating branch sign. (In section III this figure will be replaced by e^ps from (c)), (f) (f> evolution. 
Initial conditions ai t = 0, tf) — —30, Hs = 0.0148529, (p ~ 0.0699602. The following figures share these initial conditions, since 
quantum corrections are very small at this time. The evolution may be started at arbitrarily small values by evolving these 
initial conditions further backward towards the (-I-) branch vacuum. 

The resulting sources and He are plotted in Fig. |^,c,d. From the figures we see that there is no violation of NEC 
in the Einstein frame, corresponding to the fact that there is no "bounce" . This solution represents a singular collapse 
in the Einstein frame because of the linearly increasing dilaton plotted in Fig. |^f. In terms of sources this suggests 
that there is insufficient NEC violation and that the addition of conventional sources to Cm, hke radiation, which do 
not violate NEC cannot help the completion of the exit transition. 
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FIG. 2. ^Cq = — (V(7l))*, jCm = 0. See Fig. caption for details and initial conditions. 

Perhaps generically any correction violating NEC strongly enough will complete the transition. In particular, we 
consider terms modeling the form of quantum loop corrections parameterized by powers of e*^. For example, at one 
loop. 



lr4> 
2'~'q 



(2.15) 
(2.16) 
(2.17) 
(2.18) 



As long as we include only one order of loop correction, the overall coefficient of our corrections can be absorbed by 
a shift of (/), and it therefore determines the value of at which the quantum corrections begin to be important, but 
doesn't lead to qualitatively different behavior. So at this stage we choose coefficients of order unity and after these 
preliminary investigations are complete, we will choose more realistic values. 



Clearly, because of the factor of e"^ on the right hand side of equations (2.2), (2.3) this £^ can give us NEC 
violation of increasing strength as increases. We show the results of the numerical integration by presenting the 
same suite of figures as in the previous case in Fig. |^. In this case we have chosen to graph the sources for a range 
of time emphasizing the Einstein frame bounce. As hoped, the solution now proceeds into the p > region triggered 
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by increased NEC violation in the string frame. We also see the accompanying Einstein frame NEC violation and 
bounce. Notice that the bounce occurs quickly, during a few units of time, compared with the very long duration of 
the dilaton- driven phase near the (+) vacuum. 

However, inspecting the solution at late times shows the dominant terms in the equations of motion are: 

- 6i7| - Se'^f/)'* + 6iJ|(/) = (2.19) 

-6i?| + e'^<j)^ + AHl^ = (2.20) 

-6ff| + 12e'^ iJs</>3 = (2.21) 

This system has the solution, 

i/5 = f (2.22) 

approximating the nature of the true solution at late times. 

This form confirms that the solution has unbounded growth in the curvature and dilaton. Nonetheless, we have 
succeeded in our aim of completing the exit to the p > region of phase space. It may appear that we are now facing 
a new "graceful exit" problem since the equations of motion at late times are dominated by corrections, spoiling the 
expected stability of a (— ) branch. We suggest that the source of this instability is the continued NEC violation and 
will search for a cure in the next section. 

To further explore the sufficiency of NEC violation for exit we examine another generic form of one loop quantum 
correction, 

Cf = i?2 (2.24) 
pf = -imHlHs + 18ij| - 36HsHs (2.25) 
pf = IOSHIHs + bAHl + 72HsHs + 1'2hP (2.26) 
A^cf = (2.27) 

We immediately notice that for evolution strictly in the fixed point, where Hs — constant, this correction will be zero 
indicating it will not effectively contribute to NEC violation. Furthermore, this introduces higher derivatives into the 
equations of motion, which in turn introduces dangerous pathologies into the solutions |p^ . 

These pathologies come from the extra degrees of freedom in the system coming from the extra initial conditions 
that need to be imposed. These extra degrees of freedom are associated with unstable modes of the solutions which 
we regard as physically spurious, since they are solutions in which the "correction" becomes much larger than the 
terms to which it is supposed to be a small correction. Reference suggests two remedies to this problem. The first 
is "reduction of order" , in which we differentiate the large terms in the equations of motion (in this case the tree level 
terms) and use them to rewrite the higher derivatives in terms of lower derivatives, justified by the assumption the 
corrections will remain small. The second is simply to carefully chose initial conditions to avoid the unstable modes. 

Since our purpose is to make a qualitative survey of the effect of corrections and the reduction of order prescription 
is computationally prohibitive, we take the second approach and choose initial conditions so that the evolution is not 
dominated by the corrections for a reasonable span of time. Exploring these solutions confirms that corrections 
themselves do not help the exit process. The other one loop curvature squared terms, and R^'^'^^ Rf^uap yield 

sources which are constant multiples of those of R^, and as is well known, the Gauss-Bonnet combination, which does 
not contribute higher derivatives, vanishes in the equations of motion. 

But putting curvature squared terms together with corrections that do not vanish in the fixed point, for example 

Cq — 1— , does yield qualitatively different solutions which do exit, in this case a solution approaching string 

frame deSitter (with constant Hs), though again with still growing dilaton as illustrated in Fig. |[ We set our initial 
conditions near a later phase of the (+) branch vacuum to avoid instabilities. 
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Finally, we present a two loop Gauss-Bonnet correction which we choose because it represents the influence of 



curvature squared terms but does not contribute higher derivative terms. 

ijrfoB ^ g0^2^^ (2.28) 

pRls ^ e'^{-2AHl4>) (2.29) 

pRls ^ e'^(16i?|(/) + IGHsHs^ + SH^^^ + 8Hl4>) (2.30) 

A^rf^^ = e'^(24i?| + 2AHlHs) (2.31) 



Again we find that introducing this correction with an appropriate sign can complete the transition to p > as in 
Fig. m but with increasing domination by the correction terms leading to a singularity soon after the transition. We 
have also tried other combinations of sources with different coefficients and found that many of them yield solutions 
that are similar to the ones we presented, leading us to believe that our result are quite general and do not depend 
in a strong way on particular initial conditions or coefficients. 
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(e) (f) 

FIG. 4. \Cq '^^ — e'''RQg, Cm. ~ 0. See Fig. M caption for details and initial conditions. 

In summary, we have seen that generic forms of quantum corrections can complete the exit from the fixed point 
of 0] to the region of p > 0, showing the NEC violation is not only necessary, but is in some sense sufficient. The 
resulting solutions are quite varied, but we have also noted they have unbounded growth of the dilaton at late times, 
continue to be dominated by corrections and continue to violate NEC. In spite of being (— ) branch solutions they 
are still unstable. However, inducing the exit seems to be a generic property of NEC violation and not of the specific 
form of the corrections. In the next section we will attack this new exit problem, the exit from the epoch of correction 
domination. 



III. A MODEL FOR TRANSITION TO DECELERATED AND STABLE EVOLUTION 

We have seen that by using plausible forms for quantum corrections we can induce NEC violation and push the 
evolutions into a region we would like to call a completed exit. However, we have also seen that these evolutions are 
dominated by the corrections and have singular behavior unlike the desired (— ) branch solutions. We associate this 
behavior with two overlapping sources. First is the continued NEC violation itself, which tends to feed accelerated 
evolution in the Einstein frame. Second, and more directly, it is the continuing growth of that supports the strength 
of the quantum correction terms through the powers of e*^ occurring in the equations of motion. 
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Since these solutions do not spontaneously suppress the corrections we might hope that simply controlling the 
growth of the dilaton would tame the solutions. This can be don e by modeling radiation production, which can slow 
the change in the dilaton through the 3i?s(^ "friction" term in (2.4), or more directly by capturing the dilaton in a 
potential minimum. 
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-0.02(/)^e'*. Quantities plotted and initial conditions are as in Fig. [^a. Also shown are the singularity 
The right curve is for (p —> —oo and the left curve is for the value of (f> at which the solution hits 



Concentrating on the simplest of the NEC violating quantum corrections, ( 2.15| ), we are surprised to find obstacles 
to this program. Attempting to moderate the evolution through the use of a potential strong enough to affect the 
solution actually drives the solution into cusp singularities (singularities in (j) and Hs). The sources for a potential 
are: 



= -V{cl>) 



Pn 



V 



A^Cl = V'i<t>) 



(3.1) 
(3.2) 
(3.3) 
(3.4) 



We illustrate a particular case in Fig. |, with V{(j)) = - 
singularity curves. The curve on the right is the location 
their location when the solution collides with them. Adding 
redshifted away by the growing Hs. 

The source of these singularities is easy to understand 
containing the highest derivatives Hs and (j) to the terms 
we solve the equations of motion for Hs and (j) in terms of 
general, zeros in the {x,y) plane, where x = and y = 
the singularity curves: 



O.Olip'^e'^. We have added graphs of the locations of the 
of the singularities at (/) ^ — oo. The curve on the left is 
; radiation can produce similar singularities or it is quickly 

mathematically. The sources introduce additional terms 
coming from the lowest order terms in the action. When 
the lower derivatives we find a denominator which has, in 
Hs. In this particular case this leads to the equation for 



-9x^-36 e-^^*) x^ +2xy + 27x^y + 108 e'<'^^^ x^y + 2y^ + y^ = 



(3.5) 



The terms containing the e"^ are coming from the quantum correction, so that the location of the singularity curves 
is now a function of (f). 

Perversely, the singularity curves follow the solution towards small 4> and large Hs. So any attempt to tinker with 
the solution causes it to collide with them. Modeling the production of radiation produces a similar effect. While we 
do not have a direct physical interpretation of these singularities we do regard them as an indication of the general 
instability of NEC violating solutions. 
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FIG. 6. Cq = Cq^'^\ f{4') ~ ^6^(1 — tanh(8<^ -I- 36)). See Fig. ^ caption for details and initial conditions. All quantities are 
nonsingular, but peaks are out of scale to emphasize details. 

A direct approach to completing the exit transition is to assume that there exists some mechanism that shuts off 
the correction terms, and hence, NEC violation. A concrete way of modeling such an idea is simply replacing the 
quantum correction in the action with /((/>)£^ where f{(j)) is a positive constant for (j) < (po for some constant (po and 
then smoothly goes to zero, so /(</)) has the form of a smoothed step function. This successfully eliminates the loop 
corrections at late times so the dilaton may be easily captured by a potential or slowed by radiation production as 
shown in Fig. |^. In this suite of figures we have dropped the branch sign graph, since all of the following solutions are 
similar in that region of the evolution and have now chosen a range of time to graph the sources emphasizing the epoch 
when NEC violation ceases and the evolution becomes decelerated. We have also tinkered with the normalization 
of this term so the figure may be directly compared with the previous figures. While we have exactly the desired 
behavior, a function like f{(f>) will not appear by summing a few terms in the loop expansion. 



=-/(</>) (V0)4 (3.6) 

pfW = -3/(0)04 (3.7) 

P^^^^ - ~fW^' (3.8) 

A^4(*) = 12/(0)i/s0' + 3/'(0)0* + 12/(0)020 (3 9) 
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FIG. 7. £g = 1000£^ + 1000£q'*. The peaks of e^ios + ps) — 55 and pE + Pe — 55 have been cut off to emphasize the sign 
change, but all quantities are non-singular. See Fig. M caption for details and initial conditions. 

However, we have found that a complete suppression of the quantum corrections is not necessary. Suspecting that 
the instabihties are due to continued NEC violation we propose what is, perhaps, a simpler model. Since the NEC 
violation was induced in our models by one loop quantum corrections, higher loop terms can suppress the NEC 
violation once it is no longer needed if they have the correct sign. For example, a two loop contribution of the form. 



i£f ^ e^iV(bf (3.10) 
pI"^ = Se-^^^ (3.11) 
= e^'P' (3.12) 



can overwhelm the one loop NEC violation when cj) becomes large enough. Since now the scaling of different terms 
with respect to a shift in the dilaton (which determines the value of the string coupling at which various terms become 
important) is more complicated we will introduce these corrections with explicit large coefficients accounting for the 
expected large number of degrees of freedom contributing to the loop corrections. We expect, in string theory, large 
and approximately the same order of magnitude coefficients for certain one and two loop corrections, which, perhaps, 
may even be justified with some large N techniques. Taking Cq = Ci£^ + C2>Cg"^ we observe that since the sources 
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occur in the equations of motion with coefficients Cie'^ and €26^'^ respectively, and these terms will be important 
to the evolution when these coefficients are of order unity, then having Ci « C2 S> 1 leads to a situation where the 
one-loop terms become important at a smaller value of 4> (and therefore an earlier time) than the two-loop. Thus we 
can still have an era of NEC violation which is ended by the onset of the two-loop terms in a rather natural way. We 
have numerically solved the equations for a range of coefficients and observe a generic behaviour which we illustrate 
with the specific example in Fig. 0. 

We show the results of a sample evolution in Fig. 0. For the graphs of He and pE +Pe we choose a time range to 
emphasize the epoch after the bounce, where NEC violation ceases and the evolution becomes decelerated FRW. In 
Fig. ^ we show that with this form of corrections the behavior is mild enough that it is easy to capture the dilaton into 
a potential minimum. In these figures we emphasize the late phase where the dilaton is rolling around the potential 
minimum, but the early evolution is indistinguishable from Fig. M. 
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We would also like to show that these solutions are stable enough that the growing dilaton can be halted by 
introducing radiation, and that they can pass into a radiation dominated phase and be smoothly joined to standard 
cosmologies. We search for the simpl est c onsistent, but not necessarily realistic, method of producing radiation. Since 
the radiation conservation equation ( |2.5| ) can be derived from the other three equations of motion, simply placing an 
arbitrary source into ( ^.Sf ) is not satisfactory, and produces an inconsistent system. Practically, since we are using the 
equations containing the highest derivatives to integrate the system, this means that when the radiation source turns 
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on, we will begin violating the constraint equation (2.2). 

Instead we use the same ansatz used to model radiation production from the oscillation of the inflaton in slow-roll 
inflation models. We will produce the radiation from the dilaton kinetic energy. To do thi s w e introduce a coupling 
of the radiation to the dilaton by introducing a A^iZrad into (2.5) and the same term into (2.4). The natural form to 
use is A^£rad oc 4>, since this will ensure our radiation source is non-negative. We repeat for emphasis that we do not 
claim this is the actual way radiation is produced, especially since it violates the generic expectation that the dilaton 
will couple to the trace of the energy momentum tensor which vanishes for radiation. The true source of radiation 
will likely consist of the particles produced abundantly in the dilaton-driven phase. But in the spirit of this paper, it 
will serve to model the physics. 

The results are shown in Fig. ^. The produced radiation slows the dilaton to a halt, in turn suppressing the 
corrections and creating a radiation dominated evolution as is shown in Fig. In Fig. ^ we have plotted both the 
total pe + Pe and the contribution from the radiation alone pji +pr. The relation between radiation density in the 
string frame and Einstein frame is given by = e'^'^ps- 



0.20 



0.15 



0.10 




-20.0 



50.0 



40.0 



30.0 




(c) 



20.0 



10.0 



60.0 



62.0 



64.0 



(e) 



0.50 



0.40 



0.30 



0.20 



0.10 



0.00 



60.0 



61.0 



62.0 



63.0 



64.0 



Pe+P 



(b) 




(f) 



FIG. 9. Cq = 1000/:^-hlOOO£^'**, A^Crad = -100. See Fig. caption for details and initial conditions. We have superimposed 



PR +Pr, the source quantities for the radiation onto Fig. vid to show the onset of radiation domination 



These provide the first concrete examples of a completed graceful exit based on a classical evolution from an 
effective action. We have presented solutions interpolating between the inflationary (-I-) branch to decelerated (— ) 
branch evolution in which the dilaton can be captured by a potential or stopped by radiation production. Using our 
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analysis of the qualities and energy conditions required of sources to produce this transition we have arrived at an 
elusive destination. 



IV. SUMMARY AND CONCLUSIONS 

Graceful exit transition in string cosmology is not forbidden in principle as a variety of concrete examples show 
explicitly. We have verified the general arguments of showing that NEC violation is not just a necessary condition, 
but in some generic sense, also a sufficient condition. We have encountered yet another "exit problem" from a 
correction dominated evolution to a standard decelerated FRW evolution, which had to be overcome. We suggested 
that effective terms coming from higher loop corrections may do the job, and presented, for the first time, an effective 
Lagrangian whose equations of motion possess non-singular solutions interpolating between a (+) branch vacuum and 
ordinary radiation dominated FRW evolution with a fixed dilaton. 

The remaining questions concern whether specific string models produce coefficients of appropriate sign and size. 
We need to know if the one- loop terms do indeed violate NEC and whether a physical shut-off mechanism does operate 
in string theory. Note that the form of induced terms at one-loop is guaranteed, from general arguments, to be the 
one we used. 

Each of our non-singular string cosmology solutions provides a detailed description of the high-curvature phase in 
between the dilaton-drivcn inflationary phase and the FRW decelerated phase (the "string phase" ) . In all solutions we 
find that a phase rich in structure appears which is much more complicated than the one postulated in of constant 
(j) and Hs- We do not suggest taking the details of the evolution in the string phase too seriously, because the terms 
that determine those details have arbitrary coefficients. However, our examples could be taken as an illustration of 
what the real string phase may eventually look like. 
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APPENDIX: EFFECTIVE TERMS, EQUATIONS OF MOTION AND THEIR SOLUTIONS 

Before embarking on determining whether a definite string model has correction terms of the form required to induce 
the graceful exit, we have made this survey of likely forms of source terms. To do this required the variation of many 
forms of Lagrangian terms with many combinations of coefficients. Our interest was not in any one specific action. 
To make this feasible we have developed software to handle many of the aspects of the journey from action, through 
equations of motion and numerical integration and finally to graphical representation of the resulting dynamics in a 
useful form. 

The core of the process is the automated derivation of the equations of motion by varying the action, in a 
form suitable for numerical integration. Making the process difficult was the requirement that we handle essen- 
tially any form of correction term in the action. But making the prospect easier was our specializing on homoge- 
neous and isotropic solutions and the use of a tireless symbol manipulator, Mathematica [l^ . This enables us to 
use an almost embarrassingly blunt approach. First we construct a matrix to represent the metric, in our case, 
diag{gQo{t), gii{t, r), g22{t, r), gssit, r, 4>)). Then we construct the tensor quantities we need in the most direct possible 
way. We use the metric to compute the Christoffel symbols and proceed to the Riemann tensor and contractions 
thereof, all of which are stored in lists. From here we compute any required geometrical scalars for the action, which 
all emerge in a raw form, completely in terms of .goo •••533 and their various partial derivatives and any other fields. 

From here we construct the p and p and other quantities by directly varying with respect to the metric fields and 
other quantities. Our techniques for doing this were helped greatly by study of [p^ . Then we put in the cosmic time 
gauge choice 500 = ~1 and the components of the usual FRW metric a^tY'dQ, for gn, 522, 533^ Finally we replace the 
derivatives of a{t) with their corresponding expressions in terms of the Hubble parameter H = a /a. 

The very crudeness of this process enhances our belief in its correctness. Other consistency checks are possible, e.g. 
we make the redundant check that = T| = Tg . We also verify the conservation equation ( ^.5[ ) for these sources. 
They replicate many known results, e.g. the vanishing of sources from the one-loop Gauss-Bonnet combination, and 
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we can reproduce the numerical integration of other examples in the lite ratu re. On a case by case basis we check the 
accuracy of numerical integrations by verifying the constraint equation (2.2). 

Finally, we present a table of generalized sources, sufficient to construct all of the equations of motion used in this 
work, a generalized form of t he d i lato n kinetic term (Al), the Ricci scalar with arbitrary dilaton dependence (A2) 
and various combinations (A3, A4 and A5). fc is the sign of the spatial curvature. 

jr9(mf) ^ /(0)g((V0)2) (Al) 
2 



„((V0)^) ^ /('/')g((W) 

Pq 



2\ 



2 

^^^,((V0)^) ^ g((V0nr(0) ^ ^f^^)Hg'{iV4.r)<j> + /'(0)5'((V0)^)0^ + fi^)g'{{V^r)4> - 2/(0)<^^9"((V0)^)0 



((Vc/.)^) = -c/.^ 



2 



£f = i?/(0) (A2) 



ait) 

pq = ^ + yim^ + Vim + ■2Hr{<t>)^ + + /' 

a[t) 



A,/:f = + 6i7Y(0) + 3/'(0)i? 

a{t) 

Cf = R'ficP) (A3) 

Pf - ^™ + - 108/(^)i^^i^ + ISmn^ - - 72H^f' • 

a[t) a(t) a(t) 

36Hf\(j))H(j) - 36f{(j))HH 



a(i) a(t) a(t) a{t) 



a{ty a{ty 
24i/2/'(0)0 + 12/'((/))ij0 + 12/(0)i7(3) 

„2 18P/'(0) 72kH^f'((l)) , , mkf'{6)H ^ • , 

= TTT^ + rir^ + 72ifV'(0) + 7^ + 72H^nct>)H + 18/'(0)i/2 

a(t) a(t) a(t) 

£f = R'-R^.fW (A4) 

a(t) a[t) a[t) 

12Hf'{^)H^~12f{<f>)HH 

a(t) a{t) a(t) a(t) 



a{t) 



2kf 



a{tf 



a(t) a[t) ayt) 
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fl.-^ = R^''-^R^^^^f{cj>) (A5) 
ayt) ayt) 

a(t) a(f) a[t) a[tj 

a{t) a{t) 
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